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<N : ABSTRACT 

' We use data from the Hipparcos catalog and the Barbier-Brossat &; Figon (2000) 

. catalog of stellar radial velocities to test the hypothesis that the /? Pic planetesimal 

^ . disk was disrupted by a close stellar encounter. We trace the space motions of 21,497 

, stars and discover 18 that have passed within 5 pc of /? Pic in the past 1 Myr. (3 

I Pic's closest encounter is with the K2III star HIP 27628 (~0.6 pc), but dynamically 

"^H I the most important encounter is with the F7V star HIP 23693 (~0.9 pc). We calculate 

Qh| the velocity and eccentricity changes induced by the 18 perturbations and conclude that 

O ' they are dynamically significant if planetesimals exist in a /? Pic Oort cloud. We provide 

• a first-order estimate for the evolutionary state of a /? Pic Oort cloud and conclude that 

■ the primary role of these stellar perturbations would be to help build a comet cloud 

J> . rather than destroy a pre-existing structure. The stellar sample is ~20% complete and 

^ , motivates future work to identify less common close interactions that would significantly 

^ I modify the observed circumstellar disk. For future radial velocity study we identify 6 

■ " " ' stars in the Hipparcos catalog that may have approached /5 Pic to within 0.1 pc and 
therefore remain as candidate disk perturbers. 
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1. Introduction 

The dynamical mechanisms that dominate the formation and subsequent evolution of planetary 
systems may be broadly described as either endogenic or exogenic in origin. The formation of giant 
planets, and their apparent orbital migration (Goldreich & Tremaine 1980; Malhotra et al. 2000; 
Marcy et al. 2000) are probably the most significant endogenic source for modifying a planetary 
system over time. In our solar system, the giant planets have displaced a significant fraction of small 
bodies from their formation sites to either the Oort cloud or interstellar space. For circumstcllar 
dust disks observed around nearby main sequence stars, unseen giant planets are believed to produce 
the ubiquitous central depletions in the dust distributions (Roques et al. 1994; Pantin et al. 1997; 
Greaves et al. 1998; Liou & Zook 1999; Wyatt et al. 2000). 

However, young systems, such as pre-main sequence stars in clusters and the proplyd objects 
in Orion, experience significant exogenic forces (Larwood 1997). Close stellar encounters with 
other cluster members (Laughlin &; Adams 1998) and radiation pressure from nearby massive stars 
(Hcnney &: O'Dell 1999) remove disk material and can disrupt the vertical settling of circumstcllar 
dust and gas. During the evolution of our Solar System, the galactic tide and encounters with 
molecular clouds and passing stars decoupled comets from the planetary region, preserving them 
in the Oort cloud (Oort 1950; Fernandez 1997). 

The dust disk around the A5V star (3 Pic displays evidence for both types of dynamical 
mechanism. A planet may be responsible for clearing dust within ~20 AU radius of the star 
(Smith Sz Terrile 1984; Lagage & Pantin 1994), perturbing families of comets towards the stellar 
photosphere (Beust &; Morbidelli 2000) , vertically disrupting the disk at 50 AU radius (Burrows et 
al. 1995; Mouillet et al. 1997; Heap et al. 2000), and creating a few-hour duration, 0.06 magnitude, 
achromatic drop in /3 Pic's lightcurve (Lecavelier des Etangs et al. 1997). A recent and close stellar 
flyby perturbation may be responsible for producing the radially and vertically asymmetric disk 
structure at large radii (Kalas &; Jewitt 1995), as well as substructure in one disk midplane between 
500 and 800 AU projected radius (Kalas et al. 2000; Larwood & Kalas 2001). Since the age of /? Pic 
is between 8 Myr (Crifo et al. 1997) and 20 Myr (Barrado y Navascues et al. 1999), understanding 
the dynamics of objects surrounding (3 Pic could elucidate conditions and events that determined 
the early evolution of our Solar System. 

Here we will examine possible exogenic perturbations on (5 Pic. A key piece of observational 
evidence needed to confirm the stellar flyby hypothesis is to identify a perturbing star. Given im- 
proved trigonometric parallaxes and proper motions from Hipparcos (ESA 1997), we are motivated 
to test the stellar flyby hypothesis empirically. Section 2 presents our method for using Hipparcos 
and radial velocity data to trace the galactocentric motions of stars relative to /3 Pic. Section 3 
identifies stars that have passed within 5 pc of /3 during the past 1 Myr. In Section 4 we assess the 
significance of the stellar encounters for the evolution of a possible Oort cloud around (3 Pic. Sec- 
tion 5 demonstrates the plausibility of very close stellar encounters and presents a list of candidate 
perturbers for future radial velocity observations. 
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2. Search for (3 Pic stellar perturbers 

We search for candidate /3 Pic perturbers using data from the entire Hipparcos catalog and the 
Barbier-Brossat & Figon (2000) catalog of stellar radial velocities. The Hipparcos catalog gives the 
positions, proper motions and parallaxes of 118,218 stars from which we compute the galactocentric 
coordinates. Barbier-Brossat & Figon (2000) provide the mean radial velocities for 36,145 stars, of 
which 21,497 stars are also contained in the Hipparcos catalog. We combine the radial velocities 
with the Hipparcos proper motion data (Johnson &; Soderblom 1987) to trace the space trajectories 
of the 21,497 stars during the past 1 Myr. 

The observed disk substructure should be short-lived due to orbital phase mixing. The dy- 
namical models of Kalas et al. (2000) and Larwood & Kalas (2001) indicate that the perturbation 
occurred as recently as 10^ yr ago. However, we choose to extend the search 1 Myr into the past 
to allow for the uncertainties introduced by the assumptions inherent in the models. 

Key factors limiting the completeness of the search are the sensitivity constraints imposed by 
the two catalogs. In general, we expect that a significant fraction of late type stars near /3 Pic {d = 
19.3 pc) is undetected. For spectral types later than G5, Volume 1 of the Hipparcos catalog states 
a completeness limit: rriy < 7.3 -|- l.l|sin(6)|, where b is the galactic latitude. For (3 Pic, b ~ —30 
degrees, giving a completeness limit < 7.85 mag. At d = 20 pc this translates to My < 6.3 mag, 
which corresponds to the absolute magnitude of a K2V star (Cox 1999). At d = 30 pc. My < 5.5 
mag, which is the absolute magnitude of a KOV star (Cox 1999). Thus, the Hipparcos catalog 
includes a significant fraction of the A - G spectral types in a 10 pc radius volume around (3 Pic, 
but will miss K and M dwarfs [in Kalas et al. (2000) the (3 Pic disk perturber has a mass consistent 
with an M dwarf]. 

Estimating the completeness of the Barbier-Brossat & Figon (2000) catalog is considerably 
more difficult because it is not an all-sky survey, and the radial velocity information derives from 
different sources. Even though the majority of stars are < 10 mag, comparable to the sensitivity 
limit of the Hipparcos catalog, the Hipparcos catalog contains more than three times as many 
stars. The incompleteness of our sample therefore exceeds 50%, and in Section 3 we estimate an 
incompleteness of ~80%. Future missions such as FAME, SIM and GAIA will have the sensitivity 
to significantly improve the perturber search if radial velocity catalogs are also expanded. 

For P Pic's radial velocity we choose the recent measurement given by Grenier et al. (1999) 
because the published error, ajiv, is smaller than that determined by Barbier-Brossat Sz Figon 
(2000). Experimental trials using other published radial velocities for (3 Pic (e.g. Lagrange et al. 
1995) show negligible differences from the results presented below. 

The relatively short timescale considered here permits a first order estimate of the stellar 
trajectories using a straight line motion approximation. A correction for the two body interaction 
is only significant for very close (< 10^^ pc) and low relative velocity (< 10 km s^-*^) encounters. 
For each of the 21,497 stars we calculate a trajectory backward in time and determine the closest 
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approach distance, Dca (pc), to (5 Pic and the time of closest approach, tea (kyr), where tea = is 
the present. To test the vaUdity of the hnear approximation for deriving trajectories we use a fourth 
order Runge-Kutta scheme to solve the equation of motion in the Galactic potential for six of the 18 
stars that approach within 5 pc of /3 Pic. We adopt the axisymmetric Milky Way mass distribution 
model provided by Allen Sz Santillan (1991). This model^ has proven to be well adapted to the 
calculations of orbits (de Boer et al. 1997) and has been already applied to the motions of globular 
clusters and nearby stars (Allen & Santillan 1993; Schuster & Allen 1997). In the galactocentric 
cylindrical coordinates (R, 9, z) the small value of z/ R in the range of time considered here permits 
the decoupling of motions in the galactic plane and on the z axis. We check our calculations by 
reproducing the published results on the nearest approaches of stars with the Sun (Garci'a-Sanchez 
et al. 1999; Miillari k. Orlov 1996; Matthews 1994). The initial coordinates and velocities have 
been changed from the heliocentric to the galactocentric reference frame using the current lAU 
values for the LSR (0L5'i?=22O km s~^, Ri;,5ij=8.5 kpc). As expected, the agreement between the 
(Dca, tea) values derived from the straight line approximation and the ones calculated from the 
integrated Galactic orbits allows us to restrict our computations to the linear case for the selection 
of candidates. 

The main advantage of adopting the the linear approximation is that we can determine the 
influence of the errors for the input positions and velocities on the estimate of {Dca, tea) through 
a Monte-Carlo draw. We use the standard deviations provided in the Hipparcos catalog and the 
Barbier-Brossat & Figon (2000) catalog. We assume a Gaussian distributions for the errors, cen- 
tered on the average value and with dispersion equal to arms- We produce 10^ random draws that 
for each star lead to a probability distribution in the closest approach plane {Dea,tca)- 

Table 1 lists the stars that approached within 5 pc of /3 Pic during the past 1 Myr. For this sub- 
set of stars we then use SIM BAD to search the literature for recent radial velocity measurements 
that may not have been included in the Barbier-Brossat & Figon (2000) catalog. Grenier et al. 
(1999) have measured the radial velocity for HIP 19893 and give a smaller measurement error than 
that listed by Barbier-Brossat &; Figon (2000). We therefore re-calculate the space trajectory and 
errors using the Grenier et al. (1999) radial velocity, which is given in Table 1. We make the same 
correction for HIP 17378, which has a more accurate radial velocity measurement from De Medeiros 
& Mayor (1999). As with Pic, trials using different published radial velocity measurements for 
each star produce negligible changes in the final values of Dea and tea- 

3. Candidate perturbers with Dca < 5 pc 

From the 21,497 stars selected from the Hipparcos catalog and Barbier-Brossat & Figon (2000) 
catalog, 18 are found to have Dea < 5 pc with respect to /3 Pic in -10^ yr < tea < 0. Table 1 lists 



^Consisting of a spherical central bulge, a Miyamoto & Nagai (1975) disk, and a massive spherical halo. 
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the values of their respective trajectory parameters. Columns 2 and 3 in Table 2 give the maxima 
and la uncertainties for the probability density distributions derived from the Monte-Carlo draws 
for {Dca,tca)- Figure 1 shows the Dca,tca values resulting from the Monte Carlo draw for a subset 
of stars with Dca < 3 pc. Figure 2 presents the results in the form of isocontours representing 
the 68.3%, 95.4%, and 99% confidence levels for finding a star within a given region in the close 
encounter plane {Dca, tea)- An encounter probability is determined by dividing the contour area 
below a specified closest approach limit, by the total contour area in the Dca, tea plane. The 
probabilities that each star approached <1 pc and < 0.5 pc with confidence levels of 95.4% and 
68.3% are listed in columns 5 and 6 in Table 2. 

Three stars - HIP 23693, HIP 27628, and HIP 29958 - have Dca < 1 pc However, the greater 
uncertainties in the observables for HIP 29958 (Table 1) mean that there is only a 10% probability 
at the 95.4% confidence level that its closest approach was <1 pc. For HIP 23693 and HIP 27628, 
the probabilities that Dca < 1 pc are >50%. HIP 93506 and HIP 116250 are notable because 
2-3x10^ yrs before their closest approach with /3 Pic (Table 2), they passed <3 pc from the Sun 
(Garcia-Sanchez et al. 1999). 

The incompleteness of our sample is evident by comparing the number of candidate perturbers 
found by this experiment to the number expected. If we assume that (3 Pic's heliocentric distance, 
d = 19.3 pc, places it within the Solar neighborhood, then the stellar encounter frequency should 
be roughly equal to that of the Sun: N = 12 A Myr~^ for Dca <1 pc (Garcfa-Sanchez et al. 1999). 
Thus, our search produces a factor of ~5 too few perturbers. Using the Hipparcos catalog, Garci'a- 
Sanchez et al. (1999) give their empirical finding that A'" > 3.5 D^„^^ Myr~^, and they conclude 
that their sample is incomplete by at least 50%. For our Dca < 5 pc cut-off, the empirical relation 
given by Garcia-Sanchez et al. (1999) yields N > 106 Myr^^ but for /? Pic we find iV = 18 Myr-^ 
(Table 1). Again, the discrepancy indicates that onr sample is '^20% complete. Thus we expect 
that a volume- limited search realized with future stellar catalogs will reveal at least ~10^ stars 
approaching /3 Pic within 5 pc. 

To determine which candidate perturber had the greatest dynamical impact on /3 Pic, we factor 
in the relative stellar masses and velocities. We calculate the velocity impulse Av due to each stellar 
passage on both a /? Pic disk particle (r=10^ AU) and on a hypothetical Oort cloud object (r=10^ 
AU). Under the impulse approximation, the change of velocity Av of a comet relative to (3 Pic due 
to the influence of a passing star can be approximated as. 



where is the mass of the passing star, Al^a is the relative velocity, and Dca and r are the 
distances from /? Pic to the passing star and to a comet, respectively. /? denotes the angle between 
r and Dca at the time of closest approach (Oort 1950; Fernandez &; Ip 1991). 

From Av we can also estimate the change in eccentricity, Ae, using the following result from 
numerical simulations (Brunini 8z Fernandez 1996): Ae ~ 2Av/v, where the comet's orbital velocity 




(1) 
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V is v'^ = QMppic/r. Integrating Eq. 1 over /3 we compute the average perturbation of a given 
close encounter on a shell of comets at distance r. The resulting and Ae for both the maximum 
perturbation (/? = 0) and the average one at r = 10^ AU are compiled in Table 3, together with 
the 1 (J uncertainties deduced from our Monte-Carlo draws, and plotted in Figures 3a, 3b, 4a, and 
4b. For r = 10^ AU, the values for Aw and Ae are smaller by a factor determined from Eqn. 1. 
For example, with the HIP 23693 encounter, /^Vmax is a factor of 4.4 x 10^ smaller, and Ae^aa; is 
4.4 X 10^ smaller. 

Among the initial best candidates, HIP 27628 (#8) produces the greatest /^Vmax on particles 
located between itself and (3 Pic {^Vmax — 0.7 m s~^) at the time of closest approach (Table 3). 
However, its effect on the whole cloud is diminished due to its small mass and large relative velocity 
with respect to (3 Pic. Similarly, HIP 29958 (#10), with ~ O.SMq and AVca =102.1 km s"\ 
does not induce on average any notably large perturbation over the orbiting comets, even at r = 10^ 
AU. HIP 23693 (#6) remains then as the most significant perturber as it combines the smallest Dca 
in our sample, a rather small relative velocity to (3 Pic (22ib3 km s~^) and a large mass (~ I.4M0). 



4. Dynamical influence on the evolution of a /3 Pic Oort cloud 

The candidate stellar perturbers detected by our search cannot account for the ~0.005 pc close 
encounter simulated by Kalas et al. (2000) to explain the observed dust disk morphology. However, 
the values for At; and Ae at r = 10^ AU may be large enough to impact the dynamical evolution 
of objects located beyond the detected dust disk, in a /? Pic analog to the Solar Oort cloud. For 
a circular orbit at r = 10^ AU, the escape velocity for a f3 Pic comet, 176 m s~^, is significantly 
greater than the orbital velocity, 124 m s~^. Thus, only the cumulative effect of perturbers will 
dynamically modify a (3 Pic Oort cloud. On the other hand, the values for Ae are large enough 
that objects already on eccentric orbits, e > 0.99, could be stripped from the system after a single 
stellar perturbation. Moreover, objects with r > 10^ AU will experience a stronger perturbation 
that could lead to ejection. Below we determine how the Dca values for the candidate perturbers 
compare to the radius of a /3 Pic Oort cloud, and then discuss the possible evolutionary effects of 
the perturbations. 

The maximum size of a planctcsimal cloud gravitationally bound to (3 Pic is defined by the 
Roche surface of the star, at, set by the galactic tidal field (Tremaine 1993): 

Holmberg &; Flynn (2000) give the variance-weighted average from seven different studies for the 
galactic mass density, p = 0.11 MqPc~^. If Mf^pi^ = 1.75 Mq (Crifo et al. 1997), then at = 2.2 x 10^ 
AU, or ~1.1 pc. Thus, at least three of the candidate perturbers (Table 2) penetrate /3 Pic's Roche 
radius. 
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Other perturbers, though outside the Roche radius, pass close enough to temporarily exceed 
the gravitational influence of the galactic tidal field. The distance from [3 Pic where the gravitational 
forces of (5 Pic and the perturber on an Oort cloud comet are equal is given by (Miillari &; Orlov 



The last column of Table 3 gives Req for each of the candidate perturbers. In addition to the three 
perturbers that physically enter the Roche radius, HIP 25544 has Req < at- Thus a total of four 
stars out of the 18 in Tables 1 and 2 penetrate a possible Oort cloud around /3 Pic. 

In general, stars passing near an Oort cloud will either destroy a fraction of the cloud by 
sending comets into interstellar space or closer to the central star, or help build the cloud by 
increasing the periastron distances and hence the dynamical lifetimes of comets in the presence 
of planets. Stellar passages within an Oort cloud may induce comet showers which may briefly 
increase dust replenishment near the star (Weissman 1996). To determine which outcome results 
from the stellar perturbations identified here, we must first consider if the existence of a /? Pic Oort 
cloud is plausible, and if so estimate its evolutionary state. 

The four basic conditions for creating an Oort cloud around a star are (e.g. Fernandez 1997): 
1) the formation of planetesimals in a region influenced by planets, 2) the existence of a massive 
planet to dynamically pump the semi-major axes of planetesimals from their formation site to large 
distances, 3) exogenic perturbations that decouple the planetesimals from the planets by increasing 
the planetesimals' periastra, and 4) sufficient time (stellar age > relevant timescales). 

For the first condition, the existence of meter to kilometer sized planetesimals around (3 Pic 
is inferred from: a) the short lifetime of dust particles relative to (3 Pic's stellar age that implies 
a source of replenishment, probably from the coUisional erosion of larger, unseen parent bodies 
(Backman & Paresce 1993), and b) variable, transient absorption features modeled as the rapid 
sublimation of comet- like bodies near the photosphere (Beust & Morbidelli 2000). For the second 
condition, indirect evidence for a massive planet at ~1-10 AU radius is summarized in Section 1. 
The third condition, exogenic perturbations, are expected from passing stars, molecular clouds, 
and the galactic tidal field. 

For the fourth condition, wc compare (3 Pic's age, tfjpic, to the timescales required to produce 
an Oort cloud, as quantified by Trcmainc (1993). Two of the most important timescales are the 
diffusion timescale, t^, and the freezing timescale, tf. The diffusion timescale is the time for comet 
apastra to diffuse out to the Roche radius of a star at constant periastra (i.e. the time for a comet's 
energy to reach the escape energy from repeated interactions with a planet during periastron). 
The freezing timescale is the time for the galactic tide to increase comet periastra beyond the 
planet region, thereby freezing any further increase in the apastra. For a mature Oort cloud (i.e. 
well-populated and dynamically long-lived) around (3 Pic we require tf <td < tppic- 

li tf > td, then comet apastra will grow too quickly and achieve escape energy with one final 



1996): 




(3) 
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encounter with a planet. The timescales are therefore sensitive to planet mass, Mp, and semi-major 
axis, Up. If tffpic < tfi, then the system may be in the process of forming an Oort cloud, but the 
comets are probably still coupled to the planetary region. 

To satisfy td < t*, Tremaine (1993) gives the following expression: 



Mr, ( M, 



•■V 



Me - \Mq) I^IOV; K^Au) 
For tf <td, the following relation must be satisfied: 

Me" ■ \0.15Mqpc-^ J \1AU J 

The age-indepcndcnt constraint given by Eqn. 5 is significant because it generally implies that no 
dynamically stable Oort cloud will form if Mp > Msatum, even for a wide range of stellar masses. 
For example the Saturn-mass extrasolar planet candidates around HD 16141 and HD 46375 (Marcy, 
Butler & Vogt 2000) will not produce an extrasolar Oort cloud. 

In Figure 5 we plot the regions that satisfy tf < < tgpic for a range of Mp and a^, given 
MppijMQ = 1.75, p = 0.11 Mqpc-^ and tppic = 10^, 10*^ and 10^ yr. For /? Pic's estimated age 
~10^ yr, there is no combination of Mp and Up that is capable of producing an Oort cloud by the 
present epoch that is stable from ejection. In fact, a stable Oort cloud is possible only towards the 
end of P Pic's lifetime, ~10^ yr. If we assume a combination of Mp and ap that will eventually result 
in a /? Pic Oort cloud (e.g. Mp = 10 M®, Op = 5 AU), then at the present epoch comet apastra 
are only a few percent of their final values at ~10^ AU. The dust disk is detected in scattered light 
as far as ~ 2 x 10^ AU (Larwood & Kalas 2001), or ~1% of the Roche radius. Depending on how 
future observations refine our knowledge of a planetary system around (3 Pic, the observed dust 
disk may trace the early evolution of an extrasolar Oort cloud. 

In summary, objects orbiting P Pic near 10^ AU radius might still be coupled to P Pic's 
planetary region. The galactic tidal field is insufficient to decouple these objects from the planetary 
region before they attain escape energy. The stellar perturbations identified in our search neither 

eject a significant number of comets into interstellar space, nor do they cause comet showers near 
the star. The evolutionary role of these stellar perturbations is to build P Pic's Oort cloud by 
decoupling comet periastra from the planetary region. 



5. Candidate perturbers for future radial velocity study 

The most significant constraint on the completeness of this search is the lack of radial velocity 
information for nearly 100,000 stars in the Hipparcos catalog. The completion of an all-sky radial 
velocity survey, such as one proposed for the GAIA mission, has the potential to dramatically 
increase our knowledge of galactic space motions. 
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However, given only the position, distance and proper motion data from the Hipparcos catalog, 
it is possible to exclude the stars that would not approach P Pic for any physically meaningful value 
of radial velocity, and identify a sample of stars that could plausibly encounter (3 Pic. We adopt an 
iterative approach, first assuming -120 km s-i < Rv < 120 km s^^ ARv=W km s-^ to calculate 
Df-a and tea during the past 2 Myr. We select the stars with D^a < 5 pc, and repeat the procedure 
using ARv= 1 km s~^, and finally ARv= 0.1 km s~^ for the stars with Dca < 1 pc We fix the 
initial positions, distances, and proper motions to the mean values given by the Hipparcos catalog, 
rather than utilize the Monte Carlo method discussed in Section 2. 

The results are displayed in Tables 4 and 5. Tables 4 gives the maximum number of perturbers 
that can be "produced" by adjusting the radial velocity in order to get the closest possible approach 
to P Pic. A maximum of 93 stars approach /? Pic with Dca < 1 pc and tea > -1 Myr. As expected, 
very close and very recent crossing events are rare, and the more the time and distance constraints 
are relaxed, the more events are produced. The values given in Tables 4 are upper limits to the 
number of encounters experienced by P Pic, exceeding the expected rate of <1 pc encounters per 
Myr by a factor of 7.5 (Section 3). 

Table 5 gives the 22 perturber candidates with Dca < 0.5 pc during the past 1 Myr. The 
columns designate the: (1) name of the star, (2) distance of the closest possible approach, (3) 
corresponding time, (4) radial velocity that produces Dca, (5) interval of radial velocities such that 
the distance of closest approach remains below 0.5 pc, and (6) relative velocity at the closest possible 
approach. Six candidate perturbers (marked with an asterisk) have Dca < 0.1 pc and are capable 
of perturbing the observed dust disk. Table 5 both gives a specific list of objects that require 
follow-up observations to determine their true radial velocities, and demonstrates the plausibility 
of very close stellar encounters with /? Pic. 

6. Summary 

We have tested the hypothesis that P Pic experienced a close stellar fiyby by tracing the space 
motions of 21,497 stars 1 Myr into the past. We used the Hipparcos catalog for position, proper 
motion and parallax data, and the Barbier-Brossat & Figon (2000) catalog for radial velocity data. 
The completeness of our sample is ~20% and is most sensitive to B-G spectral types. Our findings 
are: 

1) In the past 1 Myr, 18 stars in our sample have had a closest approach to /? Pic Dca < 5 pc. 

2) Two stars, HIP 23693 and HIP 27628, have >50% probability that Dca < 1 pc The 
probability that Dca < 0.01 pc is negligible, meaning that these encounters cannot account for the 
large-scale P Pic disk asymmetries modeled by Kalas et al. (2000) and Larwood & Kalas (2001). 

3) The most significant dynamical perturbation on the P Pic system is due to the HIP 23693 
encounter ~356 kyr in the past. Averaged over a hypothetical P Pic Oort cloud 0.5 pc in radius. 
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the encounter induced a mean velocity change ~ 0.3 m s , and eccentricity change ~0.005. 

4) The Roche radius of (3 Pic set by the galactic tidal field is ~1.1 pc. Four stellar perturbers 
penetrate the Roche radius in the past 1 Myr and could dynamically influence planetesimals bound 
to /? Pic. 

5) We summarize evidence favoring the formation of a /? Pic Oort cloud. However, we find that 

(3 Pic is probably too young to have an Oort cloud that is decoupled from the planetary region. 
The stellar perturbations are significant in helping [3 Pic build its Oort cloud by pumping comet 
periastra away from dynamically unstable regions near planets. 

6) Wc identify a sample of 22 stars detected by Hipparcos that are candidate /3 Pic perturbers, 
but require future radial velocity observations to determine their closest approach distances. Six 
of these stars may have Dca < 0.1 and are potential disk perturbers. 
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Table 1: Input Parameters for the Stars with Peg < 5 pc 



Name 


Name 


Name 








(^1950 




'^'1950 




ll9 


50 


bl9 


'50 


Rv 




pm(cK) 






pm((5) 




"(■5) 


Spectral Type 


Notes 


(Hipparcos) 


(HD) 


(other) 


(mas) 


(mas) 


(h:m:s 


) 


(°:': 


") 




(deg) 


(deg) 


(km s^l) 


(km 




(mas) 


(mas) 


(mas) 


(mas) 






27321 


39060 


/3 Pic 


51.87 


0.51 


05:46:05.93 


-51:05:01. 


,8 


258.36 


-30.61 


19.83 


0.67 


4.65 


0.53 


81.96 


0.61 


A5V 




10798 


14412 


HR,683 


78.88 


0, 


72 


02; 


16 


:44, 


,06 


-26; 


10; 


;53, 


2 


214 


.45 


-70, 


,41 


7.00 


0, 


50 


-218.07 





.60 


444.51 





.48 


G5V 


High proper-motion Star 


17378 


23249 


(5 Eri 


110.58 


0, 


,88 


03: 


40 


:51, 


,03 


-09; 


55; 


:52, 


,9 


198 


.09 


-46, 


,00 


-6.00 


0, 


70 


-91.71 





.98 


742.23 





.92 


KOIV 


Variable of R,S CVn type 


19893 


27290 


^ Dor 


49.26 


0. 


,50 


04: 


14: 


:42. 


,85 


-51: 


36: 


;43. 


2 


259, 


.83 


-44. 


,78 


25.20 


1. 


,20 


100.57 


0, 


.46 


184.23 


0, 


.50 


F4III 


Pulsating variable Star 


19921 


27442 


e Ret 


54.84 


0. 


SO 


04: 


:15; 


:37. 


,02 


-59; 


:25: 


:18. 


,7 


270, 


.20 


-42. 


,56 


29.30 


0, 


,50 


-47.99 


0, 


.52 


-167.81 


0, 


.50 


K2IVa 


Star in double system 


22122 


30501 


GJ176.3 


48.90 


0. 


,64 


04: 


;44; 


:22. 


,88 


-50: 


09: 


32. 


,5 


256, 


.99 


-40. 


,37 


17.00 


2. 


,50 


-445.81 


0, 


.59 


-337.31 


0, 


.66 


KIV 


High proper-motion Star 


23693 


33262 


C Dor 


85.83 


0, 


,46 


05: 


;04 


:39, 


,05 


-57; 


;32: 


;25. 


,6 


266 


.03 


-36, 


,72 


-2.10 


0, 


,40 


-31.87 





.46 


117.42 





.48 


F7V 




25544 


36435 


G,72()4.1 


51.10 


0, 


52 


05: 


27' 


:03, 


23 


-60; 


27; 


15, 


,3 


269 


.31 


-33, 


,62 


12.90 


1, 


20 


-148.89 





.56 


-92.83 





.64 


G6/G8V 




27628 


39425 


3 Col 


37.94 


0, 


,57 


05: 


:49 


:11, 


74 


-35; 


:47; 


:09, 


,6 


241 


.34 


-27, 


,06 


89.10 


0, 


80 


55.74 





.50 


404.68 





.47 


K2III 


High proper-motion Star 


29568 


43162 


HR,2225 


59.90 


0, 


75 


06: 


11 


:41, 


,23 


-23; 


:50; 


:53, 


,0 


230 


.85 


-18, 


52 


22.30 


0, 


,30 


-47.06 





.45 


110.90 





.62 


G5V 




29958 




GJ9210 


32.30 


2, 


10 


06: 


16 


:04, 


,60 


-13; 


:51; 


:09, 


6 


221 


.79 


-13, 


,51 


103.80 


2, 


,00 


66.60 


1 


.91 


341.18 


1 


.64 


MOIV 


High proper-motion Star 


31711 


48189 


HR2468 


46.15 


0. 


,64 


06: 


;37: 


:28. 


,46 


-61: 


;29: 


19. 


,5 


271 


.20 


-25. 


,20 


32.30 


0. 


,30 


-50.08 


0, 


.73 


72.69 


0, 


.66 


G1.5V 


Double or multiple star 


37504 


63295 


C Vol 


24.36 


0. 


,48 


07: 


:42; 


:26. 


,97 


-72: 


:29: 


:10. 


,6 


284, 


.57 


-22. 


16 


48.10 


0. 


,50 


33.38 


0, 


.44 


15.29 


0, 


.40 


KOIII 


Star in double system 


38908 


65907 


GJ294A 


61.76 


0. 


,51 


07: 


;56; 


:51. 


,58 


-60: 


10: 


06. 


1 


273, 


.38 


-15. 


,68 


14.10 


1. 


,20 


518.38 


0, 


.63 


118.86 


0, 


.51 


GOV 


High proper-motion Star 


83990 


154577 


GJ656 


73.07 


0, 


,91 


17: 


;05 


:40, 


,03 


-60; 


:40; 


27, 


1 


329 


.82 


-12, 


22 


8.80 


1, 


20 


71.14 





.76 


589.82 





.62 


K2V 


High proper-motion Star 


89042 


165499 


L Pav 


56.32 


0, 


,68 


18: 


:05 


:47, 


,26 


-62; 


:00; 


55, 


1 


332 


.21 


-19, 


16 


29.70 


1, 


20 


-77.60 





.59 


234.68 





.44 


GOV 


High proper-motion Star 


93506 


176687 


C Sag 


36.61 


1. 


,37 


18: 


;59: 


:25. 


,89 


-29: 


;57: 


12. 


,3 


6.84 


-15. 


,35 


24.10 


3. 


,30 


-14.10 


2, 


.62 


3.66 


2, 


.00 


A2III 


Double or multiple star 


110719 


212168B 




60.24 


20.25 


22: 


:21 


:39. 


,63 


-75: 


:16: 


:11. 


,0 


314.27 


-38. 


,74 


13.00 


2, 


00 


15.18 


38.25 


61.32 


26.40 


? 


Star in double system 


114996 


219571 


7 Tuc 


45.40 


0.61 


23: 


:14; 


:31. 


,51 


-58: 


30: 


:36. 


,7 


324.30 


-54. 


,82 


18.40 


0. 


,50 


-34.93 


0.44 


79.59 


0.50 


FlIII 




116250 


221420 


GJ4340 


31.49 


0. 


,49 


23: 


30; 


:08. 


,40 


-77: 


39; 


;42. 


2 


308.32 


-38. 


,91 


26.00 


1. 


,20 


15.66 


0.51 


1.79 


0.43 


G2V 
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Table 2: Characteristics of the Closest Approach for the Stars with Peg < 5pc 



# Name tea ^co P95.4%(< Ipc) P95.4%(< 0.5pc) 

(Hipparcos) (kyr) (pc) (P68.3%(< Ipc)) (%) (P68.3%(< 0.5pc)) (%) 



1 


HIP 


10798 


-318.2t}^:^ 


1 


OO+0.51 
°°-0.40 


2 


49 


(0.00) 


0.00 


(0.00) 


2 


HIP 


17378 


-294.7l}^| 


3 


06+0-42 





00 


(0.00) 


0.00 


(0.00) 


3 


HIP 


19893 


o-i 1+35.2 
-o±.±_3g 3 


4 


04+0.19 
y^-0.14 





00 


(0.00) 


0.00 


(0.00) 


4 


HIP 


19921 


-118.9+^1/ 


3 


74+0.31 
'^-0.33 





00 


(0.00) 


0.00 


(0.00) 


5 


HIP 


22122 


-34.0li| 


2 


7^+0.56 
"J-0.18 





00 


(0.00) 


0.00 


(0.00) 


6 


HIP 


23693 


-356.3li:^ 





Q9+0.13 
^^-0.12 


54 


16 


(53.42) 


0.00 


(0.00) 


7 


HIP 


25544 


-116.0^1? 


1 


40+0.39 
^^-0.30 





00 


(0.00) 


0.00 


(0.00) 


8 


HIP 


27628 


-107.2iit? 





f:o+0.51 
^°-0.11 


70 


28 


(71.85) 


3.10 


(2.75) 


9 


HIP 


29568 


-693.2111^ 


2 


Qr+0.90 
^^-0.32 





00 


(0.00) 


0.00 


(0.00) 


10 


HIP 


29958 


-198.0li:i 


1 


QQ+2.36 
^^-0.40 


10 


21 


(12.76) 


0.70 


(0.00) 


11 


HIP 


31711 


-189.31^,^1 


3 


07+0.37 
^' -0.31 





00 


(0.00) 


0.00 


(0.00) 


12 


HIP 


37504 


-643.4li;[ 


4 


t:q+1.52 
"-•^-0.91 





00 


(0.00) 


0.00 


(0.00) 


13 


HIP 


38908 


-139.711? 


1 


Q7+0.69 
^' -0.42 





00 


(0.00) 


0.00 


(0.00) 


14 


HIP 


83990 


-303.51}!;^ 


3 


^^-1.16 





00 


(0.00) 


0.00 


(0.00) 


15 


HIP 


89042 


-397.3i^^j 


2 


20+1.21 
^^-0.21 





00 


(0.00) 


0.00 


(0.00) 


16 


HIP 


93506 


yiu.u_^28.l 


3 


rn+2.31 
^•^-1.21 





00 


(0.00) 


0.00 


(0.00) 


17 


HIP 114996 


-596.5t^|^ 


3 


04+1.35 





00 


(0.00) 


0.00 


(0.00) 


18 


HIP 116250 


-883.61^1-5 


2 


70+2.39 
'^-1.92 


3 


66 


(2.74) 


0.08 


(0.00) 
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Table 3: Dynamical Influence of the Close Encounters on Objects Orbiting (3 Pic at 10^ AU 



# 


Name 
(Hipparcos) 


Spectral 
Type 


Mass 
(Mq) 




Oca 

(pc) 


A\4o 
(km s-l) 


(10-3 m s-l) 


(10-3 m s-i) 


(10-3) 


(10-3) 




(pc) 


1 


HIP 


10798 


G5V 


0.9 


1, 


QQ+0.61 

•°**-0.40 


42.8±1.3 


33.07±30.ol 


90.52±69.02 


0.54±0.45 


1.47±0.92 


1. 


rvQ + 0.30 

^^-0.23 


2 


HIP 


17378 


KOIV 


0.8 


3, 


nc+0.42 


47.9±1.5 


4.55±0.96 


5.19±1.17 


0.07±0.02 


0.08±0.02 


2, 


qc+0.25 
■•^''-0.24 


q 
o 


HIP 


19893 




1 

1.0 


A 


04+0.19 
■^*-0.14 


1 q A-t-1 4 


lc5.oDIlIZ.iy 


ZU.OOI1IZ.40 


n qi -t-o OA 


U.043IU.U4 


9 

Z. 


cc+O.lO 
■^°-0.07 


A 
ft 


HIP 


19921 


ivzi va 


n 7 

U. I 


q 


74+0.31 
■'^-0.33 


XO.DIEU.D 


7 1 q-1-1 1 
/ . lomi. 10 


CS. 1DI1I1.41 


n 1 9-t-n (19 
u. izmu.uz 


n 1 q-i-n n9 
u. lomu.uz 


9 

Z. 


90+0.19 
■^''-0.20 





HIP 


22122 


lv± V 




9 


7R+0.56 
■'0-0.18 


Dz. / ml .0 


R ns+1 97 


7 9R-t-i Rq 


n 1 n-t-n 09 
u. lumu.uz 


n 1 9+n nq 
u. izmu.uo 


1 

1, 


R4+0.33 
■04-0.11 


6 


HIP 


23693 


F7V 


1.4 


0, 


no + 0.13 
"'^-0.12 


21.6±1.5 


308.70±97.73 


621.14+344.07 


5.03+1.59 


10.12i:5.60 


0, 


,10+0.07 
■48Io.06 


7 


HIP 


25544 




n o 
u.y 


i 


.Q+0.39 
•49-0.30 


zo.oitu.y 


c;0_l_9/1 CR 
OD.oyitz4.oD 


oi.yozt+z.oD 


n o9_i_n /in 
U.yzitu.4U 


1 qq_i_n 7n 
l.OvjitU. / U 


U. 


of.+0.23 
■°O-0.17 


8 


HIP 


27628 


K2III 


0.7 


0, 


■t>»-0.11 


83.9d=2.1 


70.82±47.74 


719.08+403.24 


1.15±0.78 


11.71±7.04 


0, 


qc:+0.31 
"^''-0.07 


9 


HIP 


29568 


G5V 


0.9 


2 


Qt-+0.90 
■yO-0.32 


11.6±0.5 


30.09±7.77 


35.36±9.88 


0.49±0.13 


0.58±0.16 


1, 


.71 +0.52 
■7II0.19 


10 


HIP 


29958 


MOIV 


0.5 


1 


•^"-0.40 


102.1±5.6 


7.42±5.27 


14.98±11.52 


0.12±0.10 


0.24±0.17 


0, 


■°"^-0.26 


11 


HIP 


31711 


G1.5V 


1.0 


3, 


Q7+0.37 
-0.31 


14.6±1.4 


18.37±3.26 


20.93±3.90 


0.30±0.05 


0.34±0.06 


2. 


9r;+0.21 


12 


HIP 


37504 


KOIII 


0.8 


4, 


"^^-0.91 


37.6±1.5 


4.39±2.71 


4.99±3.48 


0.07±0.04 


0.08±0.06 


2, 


79+0.90 
■'''-0.64 


13 


HIP 


38908 


GOV 


1.1 


1, 


Q7+0.69 
-0.42 


49.5±0.8 


22.07±8.03 


29.09±12.35 


0.36±0.13 


0.47±0.20 


1, 


nq+0-38 

■"''-0.23 


14 


HIP 


83990 


K2V 


0.7 


3 


nn+llS 

•^"'-i.ie 


60.4±1.5 


3.68±2.49 


4.33±3.72 


0.06±0.04 


0.07±0.06 


2, 


qQ + 0.72 

■"^="-0.71 


15 


HIP 


89042 


GOV 


1.1 


2, 


29+1.21 
•^^-0.21 


50.6±1.6 


12.14±4.44 


14.84±5.90 


0.20±0.07 


0.24±0.10 


1. 


97+0.67 
■^'-0.12 


16 


HIP 


93506 


A2III 


2.5 


3 


.•JU_ J 21 


39.5±5.7 


16.62±12.66 


19.37±16.94 


0.27±0.21 


0.32±0.28 


1, 


t-o+1.04 
■°°-0.55 


17 


HIP 114996 


FlHI 


1.5 


3, 


04 + 1.35 


29.3+0.8 


25.42±13.03 


49.78±28.11 


0.41±0.23 


0.81±0.47 


2, 


nq+0.70 
■"•^-0.91 


18 


HIP 116250 


G2V 


1.0 


2, 


70+2.39 
•'='-1.92 


22.9±1.7 


33.64±23.65 


55.80±41.51 


0.55±0.37 


0.91±0.73 


1, 


50+1.35 
■^°-1.09 



Table 4: Maximum number of possible perturber candidates from the Hipparcos Catalog 



Dca (pc) : 


< 0.2 


< 0.4 


< 0.6 


< 0.8 


< 1 


< 1.2 


< 1.4 


< 1.6 


< 1.8 


< 2 


t(Dca)>-100 
































t(Dca)>-200 

















1 


1 


2 


2 


2 


t(Dca)>-300 














2 


3 


3 


4 


4 


5 


t(Dca)>-400 








1 


1 


5 


7 


8 


10 


10 


11 


t(Dea)>-500 


1 


2 


3 


6 


10 


14 


20 


22 


24 


25 


t(Dea)>-600 


3 


5 


11 


16 


24 


28 


39 


41 


44 


45 


t(Dea)>-700 


5 


12 


21 


29 


43 


49 


64 


68 


81 


84 


t(Dca)>-800 


12 


23 


33 


45 


64 


73 


91 


96 


111 


119 


t(Dca)>-900 


14 


27 


40 


56 


77 


90 


110 


117 


138 


148 


t(Dc„)>-1000 


20 


36 


51 


70 


93 


111 


133 


141 


164 


174 
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Table 5: Perturber candidates for future radial velocity studies 



Name 




tea 






<0.5 pc) 


AVea 


(Hipparcos) 


(pc) 


(kyr) 


(km s~^) 


(km 




(km s~^) 


HIP 1837* 


0.025 


-537.0 


30.6 


[ 29.5 


31.7] 


37.7 


HIP 6343 


0.325 


-959.0 


89.9 


[ 87.8 


92.1] 


78.8 


HIP 6913 


0.229 


-560.0 


88.2 


[ 84.8 


91.9] 


100.9 


HIP 11324 


0.127 


-998.0 


39.7 


[ 38.9 


40.7] 


28.9 


HIP 11558 


0.374 


-680.0 


40.8 


[ 40.2 


41.5] 


41.6 


HIP 12960* 


0.007 


-732.0 


103.7 


[ 99.1, 


108.7] 


90.6 


HIP 14007 


0.327 


-910.0 


50.3 


[ 49.1 


51.5] 


34.9 


HIP 14116* 


0.064 


-799.0 


91.1 


[ 87.1 


95.3] 


74.7 


HIP 22817* 


0.082 


-725.0 


98.8 


[ 96.3, 


101.2] 


97.5 


HIP 26067 


0.490 


-283.0 


110.1 


[ 105.0 


116.1] 


90.3 


HIP 30260 


0.480 


-908.0 


58.1 


[ 56.8 


59.4] 


37.6 


HIP 33973* 


0.077 


-696.0 


80.2 


[ 77.8 


82.7] 


63.3 


HIP 42215 


0.154 


-451.0 


78.7 


[ 75.9 
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Fig. 1. — Positions in the closest approach plane {Dca,tca) from the random Monte-Carlo draw. 
From the 18 candidate perturbers, we show a subsample of 8 stars that have D^a < 3 pc. Each color 
corresponds to a different star, which have number labels corresponding to column 1 in Table 2 
and Table 3. Crosses mark the maxima of the probability distributions for {Dca,tca)- The spread in 
the distribution of points reflects the initial uncertainties on the stellar proper motions, parallaxes, 
and radial velocities. The horizontal dashed line marks -D = 1.0 pc (~ /3 Pic's Roche radius) and 
the dotted line marks D = 0.5 pc. 
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Fig. 2. — Figure 2b on next page. 
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Fig. 2. — Isocontours for the Monte-Carlo distributions in the closest approach plane {Dca,tca) 
for the 18 stars with Dca < 5 pc. Bold contour: 68.3% confidence level, thin contour: 95.4%, 
and dashed contour: 99% confidence level. The contours outline the decreasing probability of 
finding the star at a more specific location at a given time due to uncertainties in proper motion, 
radial velocity, and parallax. Each cross marks the position of the maximum of the probability 
density distribution in the Monte Carlo simulation (the size of the crosses has no significance). The 
probabilities that are given in columns 5 and 6 of Table 2 are calculated by dividing the area for 
each contour below the 1 pc and 0.5 pc lines by the total area enclosed by the contour. Figure 2b 
is an enlarged version of Figure 2a for the four candidate perturbers with the smallest Dca- 
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Fig. 3. — Velocity changes induced by the candidate perturbers on hypothetical objects orbiting 
P Pic with semi-major axis 10^ AU (equation 1 and data from Table 3). Fig. 3a (top) shows the 
mean velocity change for particles distributed in a shell with radius 10^ AU. Fig. 3b (bottom) 
plots the maximum velocity change induced by a given stellar encounter (i.e. for objects that lie 
closest to the perturbing star). Error bars are 1 a (68.4% confidence level). Stars are identified 
with numbering given in Tables 2 and 3. 
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Fig. 4. — Eccentricity changes induced by the candidate perturbers on a hypothetical objects 
orbiting /3 Pic with semi-major axis 10^ AU. Fig. 4a (top) shows the mean eccentricity change for 
particles distributed in a shell with radius 10^ AU. Fig. 4b (bottom) plots the maximum eccentricity 
change induced by a given stellar encounter (i.e. for objects that lie closest to the perturbing star). 
Error bars are 1 a (68.4% confidence level). Stars are identified with numbering given in Tables 2 
and 3. 
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Fig. 5. — Relationship between planet parameters and Oort cloud formation timescales. Open 
triangles plot the constraint for the diffusion timescale, t^, assuming three different stellar ages 
(Eqn. 4). The triangles point upwards to indicate the r egions of planet parameter space 
capable of diffusing comets to the (3 Pic Roche radius within the specified age (t^ < t^). Below this 
boundary, there has been insufficient time to produce a cloud of comets with large apastra (~ 10^ 
AU). The timescale for freezing the semi-major axis of a comet, tf, is shown by downward-pointing 
solid triangles (Eqn. 5). Below the solid triangles the galactic tidal field decouples the comets 
from the planet region before ejection energy is reached {tf < td). Comets are ejected if the planet 
parameters lie above the solid triangles. A mature Oort cloud can form around (5 Pic only after 
10^ yr when tf <td< tppi^. 



